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Abstract 

We discuss the basic transport phenomena in gases and plasmas obeying 
the g-nonextensive velocity distribution (power-law). Analytical expressions 
for the thermal conductivity (K q ) and viscosity (rjq) are derived by solving 
the Boltzmann equation in the relaxation-time approximation. The available 
experimental results to the ratio K q /r] q constrains the ^-parameter on the 
interval 0.74 < q < 1. In the extensive limiting case, the standard transport 
coefficients based on the local Gaussian distribution are recovered, and due 
to a surprising cancellation, the electric conductivity of a neutral plasma is 
not modified. 
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I. INTRODUCTION 



In the last few years, an increasing attention has been paid to possible nonextensive 
effects in the fields of thermodynamics and statistical mechanics. The main motivation is 
the lack of a comprehensive treatment including gravitational and Coulombian fields for 
which the assumed additivity of the entropy present in the standard approach is not valid 
[1-2]. Inspired on such problems, as well as in the traditional ensemble theory, Tsallis || 
proposed a remarkable g-parameterized nonextensive entropic expression which reduces to 
the extensive Gibbs- Jaynes-Shannon entropy in the limiting case q = 1 (see M for a regularly 
updated bibliography on this subject). 

Later on, the first attempts exploring the kinetic route associated to this nonexten- 
sive approach appeared in the literature [5-10]. The original kinetic derivation advanced by 
Maxwell [TT| was generalized to include power law distributions as required by this enlarged 
framework. In particular, it was shown that the equilibrium velocity g-distribution 

f (v) = B q [l-(l- q )^}^\ (1) 

is uniquely determined from two simple requirements ||: (i) isotropy of the velocity space, 
and (ii) a suitable nonextensive generalization of the Maxwell factorizability condition, or 
equivalently, the assumption that F(v) = f(v x )f(v y )f(v z ). The quantity B q above denotes 
the g-dependent normalization constant whose expression on the interval l/3<g<lis 

where n is the particle number density, m is the mass and T is the temperature. As expected, 
in the limiting case q — 1, such expressions reduce to the standard Maxwellian ones ||. 

More recently, the kinetic foundations of the above distribution were investigated in a 
deeper level through the generalized Boltzmann's transport equation || 

df df F df 
ot or m av 
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where C q is the g-nonextensive integral source term measuring the change in / due to 
collisions. This Boltzmannian like approach incorporated the nonextensive effects using two 
different ingredients. First, a new functional form to the kinetic local gas entropy, and, 
second, a nonfactorizable distribution function for the colliding pairs of particles whose 
physical meaning is quite clear: the Boltzmann chaos molecular hypothesis is not valid in 
this extended framework. It was also shown that the kinetic version of the Tsallis entropy 
satisfies an ifq-theorem, and more important still, the g-parameterized class of power law 
velocity distributions emerged as the unique nonextensive solution describing the equilibrium 
states. 

In this paper we go one step further by computing the basic nonextensive transport 
coefficients for dilute gases and plasmas. Due to the intricate form of the new collisional 



term, we adopt here the relaxation-time model proposed by Bhatnager, Gross and Krook [12 



because it permits an exact mathematical treatment. In spite of its limitation, the relaxation- 
time model (henceforth BGK model) has been proved quite useful from a methodological 
viewpoint because it yields the correct answer to the problem in a first approximation, 
and as such, it can guide us to the correct nonextensive results which must rigorously be 
obtained using the full g-transport equation. As we shall see, analytical expressions for 
the transport coefficients of a dilute gas are readily obtained with basis on the local q- 
nonextensive distribution. In particular, we show that the available experimental results to 
the ratio K q /i] q constrains the g-parameter on the interval 0.74 < q < 1. As expected, in 
the extensive limiting case, the classical expressions of the transport coefficients are readily 
recovered. However, due to a surprising cancellation, the standard electric conductivity of 
a dilute neutral plasma is not modified. The analysis presented here is also compared with 
independent calculations based on rather different approaches [f7yT0|. 
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II. NONEXTENSIVE TRANSPORT COEFFICIENTS 



Let us now consider a dilute gas whose particles are acted by external forces. In the 
BGK approximation the g-transport equation (|3|) reads : 

df df F df f-f 

■£ + v -jf + --7T = - ' 4 

ot or m (TV t 

where F is the external force and m is the mass of the particles. In this approximation, / 
is the local equilibrium distribution function and r is a relaxation time which is a number 
of the order of the collision time. We recall that under stationary conditions the distribu- 
tion function does not depend on time, however, the concentration, and the temperature, 
are local quantities. The whole system is assumed to be out but close to the local sta- 
tionary equilibrium state fl3| . In particular, this means that the resulting nonequilibrium 



distribution, /(r, v), is only slightly different from the equilibrium unperturbed stationary 
distribution, fo{r, v ) ; and can be approximated as 

/(r,v) = / (v,r) + </(v,r), \g\ < /„. (5) 

In what follows we first analyze the effects of the nonextensive distribution on the thermal 
conductivity K q . It will be assumed that the heat flux is macroscopically governed by the 
classical Fourier law 

q e = -K q VT, (6) 

while, microscopically, it is expressed as the usual average value of the perturbed kinetic 
energy flow ||14| . [T5| 

q e = -m J v 2 vgd 3 v. (7) 
As one may check, in the absence of external forces, the BGK model leads to 

9 = -tv • — . 8 

OT 
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Substituting (||) into (0), and rewriting the resulting expression in spherical coordinates 
(v, 9, 4>), one obtains 



] roc /*7T / 1 27T 

[, -~m / v 4 rdv / sin6d6 / v 
2 Jo J o Jo 



OT 



(9) 



Before proceed further, we note that the orthogonality relation 

4tt 



7T /*27T 

o jo 



ViVj sin 



— w Si 



(10) 



where j = x, y, z, may be used to show that the angular integral can be written as 



r2n 


dfo 


/ V 

/o 


v ' TT 
or 



4?r 2 df 
d(f) = —v — , 
3 OT 



(11) 



with equation (|S|) simplifying to 



2tt 



q e = — — m / v b T—dv. 
3 JO OT 



(12) 



As remarked earlier, the concentration n and the temperature T in the local g-distribution 
are spatially dependent while the pressure remains constant. Thus, as a consequence of the 
extended virial theorem ||, one may write n(r)T(r) = constant, and combining that with 
the spatial gradient of the unperturbed distribution, the heat flow can be rewritten as 

2vr 



Z7T r°° 
qe = -— mT B q v 6{[l_(l_ g ) 
3 Jo 



2A; B T J 



2k B T 



5 1 r / \ mv n1 /n >,9T , 

1 - (1 - o) — cfc. 

2T [ v HJ 2k B T i ot 
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Now, comparing with the Fourier law (Q), we obtain an integral expression for the nonex- 
tensive thermal conductivity 

^ 8 [l-(l-g)- 



w r 8[1 _ fl y |g/Md 

9 3k R T q Jo [ [ 2knT ^ 



q 3k B T 



9 Jo L V H> 2k B T ] 



3T jo zKg, 
The above integrals can easily be evaluated. We find 
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(14) 



T 2 



Jo 1 K HJ 2k B T l 
1 (2k B T\ 9/2 IOStt 1 ^ / 1 \ 1/2 r(^ 



2 \ m ) (9-7g)(7-5g)(5-3g)(3-g) \l-qj r(| + ^ 



(15) 



and, 



°°- 6 [l-(l-g)^] 1/(1 ^ 



1 / 2^T\ 7/2 SOttV^ / 1 x 1/2 r( _i_) 

2\ rn ) (9-7g)(7-5g)(5-3g)(3-g) Vl-J r(| + T: 1 



(16) 



As one may check, inserting these integrals and the normalization constant B q into equation 



(|l4l), we obtain 

5 nrk 2 B T 



K. 



q 2 m 



(17) 



_(9-7g)(7-5g)_ 

The nonextensive parameter q in the above expressions are restricted to positive values with 
q 7^ 9/7,7/5. Note also that expression ( jT7[) reduces to standard value in the extensive 
limiting case [ |i5| , |i6| 



K\ = h - U ^L. (18) 
2 m K ' 

It thus follows that the ratio between the nonextensive thermal conductivity and the exten- 
sive value, Kg/Ki, is dependent on the Tsallis' thermostatistics through parameter q (see 
Fig. 1). Note that for values of q smaller than unity, the corresponding coefficient K q can 
be much smaller than the standard result (nearly Ki/3 for q ~ 0.7), while for q > 1 it may 
increases without limit. 

Let us now consider the nonextensive effects on the viscosity coefficient. In this case, 
the particles of the gas do not move with the same velocity. In a simplified treatment, 
the particles have constant mean velocity u x in the x direction with the magnitude of u x 
depending only on z, that is, u x = u x {z). Thus, neglecting bulk viscosity, the Navier-Stokes 
stress reduces to 



p. = -%^f>. (is) 

A complementary hypothesis is that the collisional state tend to produce a local equilibrium 
distribution relative to the moving gas with mean velocity u x . Under such conditions, the 
stationary g-distribution becomes 

fo(v x - u x (z), vy, v.) = B q [l-(l- q)^] 1/{1 - q \ (20) 

where U x = v x — u x (z), U y = v y , U z = v z . Since the resulting distribution / is stationary 
but also depends on z (the direction of the gradient velocity), the BGK equation (f|) now 
reads 

9f f — fo / 01 x 

v 'Tz = -—■ (21) 

The perturbed distribution function obeys 

9 = - TV ^ ( 22 ) 

so that 

df du x (z) 

f = h+TVz w x —- (23) 

In order to calculate the component P zx of the stress tensor, we consider its kinetic definition 







P {j = J dUUiUjf. (24) 



Following standard lines we see that (for details see Ref. [16 



P zx = m[ d"Ur^Lu!U x ^^, (25) 
J dU x z dz K J 



and comparing this result with the macroscopic equation (|I9|) we find 

dfo 



Vq = -mr J £U^j-U 2 z U x , (26) 



where 



= -^[1 - (1 - q)—f^mU x (27) 
dU x k B T [ V y; 2A; B T J x V ; 



Replacing (£7]) into (p26[ ) it follows that 

% = mrA /// dZ7,dl/ lf dtr B [l - (1 - t)^<Pl + + U 2 z )]^U 2 z Ul (28) 
and performing the elementary integrals 

„ - ^L B Y 2 (^Y 2 (29) 

"'-k B T B '{ m ) (7-5 5 )(5-3,)(3-,)l,l- 5 J r(i + i y ^ 



Finally, substituting the expression of B q and cancelling out the common factors, the vis- 
cosity coefficient assumes the rather simple form 

2 

V q = ^_ 5 ^ nrk B T. (30) 

This expression is valid only for q 7^ 7/5, and as expected, in the limit q — > 1 it reduces to 
the standard extensive result |15[ 



771 = nr/c^T. (31) 

In Figure 2 we show the ratio r) q /rji as a function of the nonextensive parameter. This 
plot is similar to what happens with the thermal conductivity (see Fig.l), and as expected, 
only for q — 1 the standard Maxwellian result is recovered. Note also that combining 
equations flT7| ) and (|30"D we also obtain a quite simple expression to the dimensionless ratio, 



A, involving the transport coefficients, K q and T] q , namely 



mK„ 5 

A = q - = . 32 

Ma 9-7g v ; 



Such a quantity (sometimes called Eucken's ratio |17| ) plays an important experimental 



role because the unknown collision time appearing in both transport coefficients cancels out 
in this ratio. In Figure 3, we plot A as a function of q. Note that A does not depend 
on the temperature, and for q = 1 it reduces to the extensive result A = 2.5 (ff4j. As 



widely known, this pure number lies experimentally within the range 1.3 < A < 2.5 [16,18 



Therefore, the allowed values to the nonextensive parameter is restricted on the interval 
0.74 < q < 1, thereby showing that the Maxwellian result is only marginally compatible 
with such measurements. 

It is worth mentioning that the departure from standard Maxwellian prediction is usually 
taken as an indication that the energy of the molecules must include another forms than 
mere kinetic energy of translation, or equivalently, an anomalous specific heat . However, 
as we have seen, a possible alternative explanation (at least for monatomic gases) is provided 
by the existence of nonextensive effects associated to Tsallis' thermostatistics. 

For completeness, we now investigate the possible nonextensive effects on the electric 
conductivity of a dilute neutral plasma. As before, we consider the BGK approximation 
assuming that the external electric field is uniform and points to the positive z direction, 
that is, E = Ez. The relevant macroscopic equation is now the Ohm's law 

j z = a q E, (33) 

where a q is the electric conductivity. Microscopically, the electric current density is defined 
as the statistical average 1]T4|,|T5] 



: / d 3 vfv n , (34) 



where the subscript n means that the charge flow is directed along the normal to the cor- 
responding surface element. Note that only the component j z does not vanish since for an 
homogeneous and isotropic medium the current j is parallel to the electric field. As one may 
check, from equations ([!]), @ and fl34|) it follows that 



e2rEB " j dM l- { l- q) ^/(^ v l (35) 



k B T J y 2ksT 

Performing the integration and comparing to the Ohm's law the resulting electric conduc- 
tivity is 



e 2 r „ v2f2k B T\ 5/2 2 ( 1 \ 1/2 T(^) 



'q 



k B T * \ m J {3-q)(5-3q)\l-qJ r(i + ji-) ' 



or still, inserting the value of B q (see Eq. (||)) 



ne 2 r 



a, 



q 

m 



(37) 



Therefore, within the BGK approximation, the above cancellation results that the electric 
conductivity does not depend on the nonextensive g-parameter. 

III. CONCLUSION 

In this paper we discussed some possible nonextensive effects on the transport phenomena 
in gases and plasmas. It was shown that the theoretical analysis based on the relaxation- 
time approximation lead to simple analytic expressions to the heat conduction and viscosity 
coefficients, while the electric conductivity was not modified. 

It is interesting to compare our results (and the overall approach) with independent 
analyzes and some previous expressions to nonextensive transport coefficients appearing in 
the literature 0,0. We first notice that our main results are quite different from that ones 



obtained by Boghosian ]7| . Probably, the basic reason comes from the fact that in his paper 
it was advocated a different choice to the BGK operator, namely, C q = —r~ 1 (f q — f Q q ). 
Note that it depends nonlinearly on the nonextensive parameter and should be compared 
with the standard BGK operator assumed in the present paper (see our equation (|])). In 
the treatment adopted here, the nonextensive effects appear only implicitly through the 
power law equilibrium distribution. In addition, Boghosian (and also the authors of Ref. 
||10|| ) assumed that the pressure in the perturbed gas does not remain constant thereby 
giving rise to an anomalous heat conduction (T q in their notation). In this concern, our 
approach may be considered more conservative. As a matter of fact, using the extended virial 
theorem 0, we have assumed constant pressure as usually done in the BGK approximation 
so that the anomalous heat conduction coefficient, T q , is absent. Another crucial difference 
among these works is that three kinds of averaging methods has been considered. Boghosian 
and Potiguar et al. []TU[ employed unnormalized and normalized g-expectations values, 
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respectively, while we have computed the averages in the usual manner. In this way, it is not 
surprising that although obtaining coefficients K q and r\ q proportional to the temperature, 
the behavior of the dimensionless ratio A oc K q /r] q as a function of the q parameter is quite 
different (compare our equation ([32]) with equations (83) and (40) in the quoted papers). 
Since this ratio can be measured in laboratory, hopefully, these different approaches, as well 
as their basic underlying assumptions will be experimentally confronted in the near future. 
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FIGURES 




FIG. 1. Nonextensive thermal conductivity. We show the dimensionless ratio K q /K\ as a 
function of the q parameter. The normalizing quantity K\ corresponds to the Maxwellian expression 
obtained in the limit q = 1. Note that for q — > 1.3 the ratio K q /K% increases with no limit. 




FIG. 2. Nonextensive viscosity coefficient. The plot shows the dimensionless ratio rjq/rjx, where 
rji denotes the Maxwellian result, as a function of the nonextensive parameter q. For q — > 1.4, this 
ratio goes to infinity. 
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FIG. 3. The A parameter. The experimental dimensionless Eucken's ratio, A = mK q /kBTjq, is 
shown as a function of the nonextensive parameter q. The shadowed region shows the measured 
interval 1.3 < A < 2.5, for which the q parameter is restricted to 0.74 < q < 1. Note that the 
Maxwellian prediction is only marginally compatible with such experiments. 
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